Vertex labeling and routing in expanded 
Apollonian networks 

Z. Zhang, F. Cornelias, G. Fertin, A. Raspaud, L. Rong, S. Zhou 



1 



Zhongzhi Zhang 

Department of Computer Science and Engineering and 
Shanghai Key Lab of Intelligent Information Processing, 
Fudan University, Shanghai 200433, China 
zhangzz@f udan . edu . cn 

Francesc Cornelias 

Dep. de Matematica Aplicada IV, Universitat Politecnica de Catalunya 
Av. Canal Olfmpic s/n, 08860 Castelldefels, Barcelona, Catalonia, Spain 

comellas@ma4 . upc . edu 

Guillaume Fertin 

LIN A, Universite de Nantes 
2 rue de la Houssiniere, BP 92208, 44322 Nantes Cedex 3, France 
f ertinOlina.univ-nantes . fr 

Andre Raspaud 

LaBRI, Universite Bordeaux 1 
351, cours de la Liberation, 33405 Talence Cedex, France 
raspaud@labri . f r 

Lili Rong 

Institute of Systems Engineering, Dalian University of Technology 
Dalian 116024, Liaoning, China 
llrong@dlut . edu . cn 

Shuigeng Zhou 

Department of Computer Science and Engineering and 
Shanghai Key Lab of Intelligent Information Processing, 
Fudan University, Shanghai 200433, China 
sgzhou@f udan . edu . cn 



2 



Abstract 



We present a family of networks, expanded deterministic Apollo- 
nian networks, which are a generalization of the Apollonian networks 
and are simultaneously scale-free, small-world, and highly clustered. 
We introduce a labeling of their vertices that allows to determine a 
shortest path routing between any two vertices of the network based 
only on the labels. 

1 Introduction 

In these last few years there has been a growing interest in the study of 
complex networks [21 1201 EH HI] , which can help to describe many social, 
biological, and communication systems, such as co-author networks [30] . 
sexual networks [2H], metabolic networks [21], protein networks in the cell 
[26J, Internet [2H], and the World Wide Web jlj. Extensive observational 
studies show that many real-life networks have at least three important com- 
mon statistical characteristics: the degree distribution exhibits a power law 
tail with an exponent taking a value between 2 and 3 (scale-free); nodes hav- 
ing a common neighbor are far more likely to be linked to each other than 
are two nodes selected randomly (highly clustered); the expected number of 
links needed to go from one arbitrarily selected node to another one is low 
(small-world property). 

These empirical findings have lead to a new kind of network models [2J 
l2*TH EH! EE] • The research on these new models was started by the two semi- 
nal papers by Watts and Strogatz on small-world networks [32] and Barabasi 
and Albert on scale-free networks [7j. A wide variety of network models and 
mechanisms, including initial attractiveness [21], nonlinear preferential at- 
tachment [2H], aging and cost [3], competitive dynamics [Hj, edge rewiring [3] 
and removal [TH], duplication |T2], which may represent processes realistically 
taking place in real-life systems, have been proposed. 

Recently, based on the classical Apollonian packing, Andrade et al. intro- 
duced Apollonian networks [B] which were simultaneously proposed by Doye 
and Massen in [22J. Apollonian networks belong to a deterministic type of 
networks studied earlier in Refs. [HI El EH El E3 which have received much 
interest recently [101 E31 EZ1 EE] ■ Two-dimensional Apollonian networks 
are simultaneously scale-free, small-world, Euclidean, space filling, and with 
matching graphs [HI EH]- They may provide valuable insight into real-life net- 
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works; moreover, they are maximal planar graphs and this property is of par- 
ticular interest for the layout of printed circuits and related problems [UJ 0U] . 
More recently, some interesting dynamical processes [10J EH [2H], such as 
percolation j^U], epidemic spreading [3U|, synchronization [37|, and random 
walks [23] , taking place on these networks have been also investigated. 

Networks are composed of vertices (nodes) and edges (links) and are very 
often studied considering branch of discrete mathematics known as graph 
theory. One active subject in graph theory is graph labeling [21]. This is 
not only due to its theoretical importance but also because of the wide range 
of applications in many fields [TP] , such as x-rays, crystallography, coding 
theory, radar, astronomy, circuit design, and communication design. 

In this paper we present an extension of the general high dimensional 
Apollonian networks [HI 122 EH] which includes the deterministic small-world 
network introduced in [36]. We give a vertex labeling, so that queries for the 
shortest path between any two vertices can be efficiently answered thanks 
to it. Finding shortest paths in networks is a well-studied and important 
problem with also many applications [1. Our labeling may be useful in 
aspects such as network optimization and information dissemination, which 
are directly related to the problem of finding shortest paths between all pairs 
of vertices of the network. 

2 Expanded Apollonian networks 

In this section we present a network model defined in a iterative way. The 
model, which we call expanded Apollonian network (EAN), is an extension of 
the general high dimensional Apollonian network [U] [221 EH] which includes 
the deterministic small-world network introduced in 

The networks, denoted by A(d, t) after t iterations with d > 1 and t > 0, 
are constructed as follows. For t = 0, A(d, 0) is a complete graph (or 
(d + 2)-clique). For t > 1, A(d, t) is obtained from A(d, t — 1). For each of 
the existing subgraphs of A(d, t — 1) that is isomorphic to a (d+ l)-clique 
and created at step t — 1, a new vertex is created and connected to all the 
vertices of this subgraph. Figure Q shows the network growing process for 
the particular case where d = 2. 

Let n v (t) and n e (t) denote the number of vertices and edges created at 
step t, respectively. According to the network construction, one can see 
that at step t« (tj > 1) the number of newly introduced vertices and edges 
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Figure 1: Illustration of a growing network in the case of d = 2, showing the 
first two steps of growing process. 



is n v (ti) = (d + 2){d + l)* 1 " 1 and n e (t;) = (d + 2)(d + l) u . From these 
results, we can easily compute the total number of vertices N t and edges E t 
at step t, which are N t = i^M^+d^ and Et = {d + 2 )(d + 1) 2(d+ f/ rf " 2 , 
respectively. So for large t, the average degree k t = ^ is approximately 
2(d + l). 

This general model includes existing models, as listed below. 

Indeed, when d — 1, the network is the deterministic small- world net- 
work (DSWN) introduced in [SH] and further generalized in DSWN is 
an exponential network, its degree distribution P(k) is an exponential of a 
power of degree k. For a node of degree k, the exact clustering coefficient 
is |. The average clustering coefficient of DSWN is In 2, which approaches 
to a high constant value 0.6931. The average path length of DSWN grows 
logarithmically with the number of network vertices [39J. 

When d > 2, the networks are exactly the same as the high-dimensional 
Apollonian networks (HDAN) with d indicating the dimension jHJ 1221 CUB EH] ■ 
HDAN present the typical characteristics of real-life networks in nature and 
society, and their main topological properties are controlled by dimension d. 
They have a power law degree distribution with exponent 7 = 1 + ln ^ 1 ^ 
belonging to the interval between 2 and 3 jHl 1221 E3 EH] ■ For any individual 
vertex in HDAN, its clustering coefficient C(k) is also inversely proportional 
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to its degree k as C(k) = k ^ k _^ ■ The mean value C of clustering coefficient 
of all vertices in HDAN is very large and is an increasing function of d. For 
instance, in the special cases where d = 2 and d = 3, C asymptotically reaches 
values 0.8284 and 0.8852, respectively. In addition, HDAN are small worlds. 
The diameter of HDAN, defined as the longest shortest distance between all 
pairs of vertices, increases logarithmically with the number of vertices. So, 
the EAN model exhibits a transition from an exponential network (d = 1) 
to scale- free networks {d > 2). 

3 Vertex labeling 

Vertex labeling of a network is an assignment of labels to all the vertices in 
the network. In most applications, labels are nonnegative integers, though 
in general real numbers could be used In this section, we describe a way 
to label the vertices of A(d, t), for any d > 1 and t > 0, such that a routing 
by shortest paths between any two vertices of A(d, t) can be deduced from 
the labels. We note that a more general result on shortest paths routing 
of graphs with given treewidth is given in However, here we address 

the more specific case of the expanded Apollonian networks A(d, t). In what 
follows, we will denote L(v) as the label of vertex v, for any vertex v belonging 
to A(d,t). 

Here the labeling idea, inspired from is to assign to any vertex v 
created at step t > 1 a label of length t, in the form of a word of t digits, 
each digit being an integer between 1 and d + 2 (the vertices obtained at step 
t = 0, i.e. the vertices of the initial (d + 2)-clique A(d, 0), are assigned a 
special label). More precisely, the labeling of any vertex v of A(d, t) is done 
thanks to the following rules: 

• Label the vertices of the initial (d + 2)-clique A(d, 0) arbitrarily, with 
labels 1', 2'...(d+ 2)'. 

• At any step t > 1, when a new vertex v is added and joined to all 
vertices of a clique K d+1 : 

1. If v is connected to d+ 1 vertices of the initial (d + 2)-clique, then 
L{v) = I, where I' is the only vertex of the initial (d + 2)-clique 
that does not belong to this (d + l)-clique. 
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2. If not, then v is connected to Wi,w 2 ■ ■ -Wd+i, where at least one 
of the WiS is not a vertex of the initial (d + 2)-clique. Thus, any 
such vertex has a label L(w%) = si,iS2,i ••• Sk,i- W.l.o.g., let w\ 
be the vertex not belonging to the initial (d + 2)-clique with the 
longest label. In that case, we give vertex v the label L(v) defined 
as follows: L(v) = a ■ L(w\), where l<a<<i + 2is the only 
integer not appearing as first digit in the labels of W\, tu 2 • • • Wd+ii 
that is a = {1, 2 . . . d, d + 1, d + 2}/ ufij 1 (the fact that a is 
unique will be proved by Property 1 below). 

Such a labeling is illustrated in Figure 121 In the upper part of this figure, 
we label the vertices of A(d,t), for d — 1 and up to t = 3. We see that 
vertex u, created at step 1, has label L(u) = 2 because it is not connected to 
vertex 2' of the initial 3-clique (triangle). Vertex w is is not connected to any 
vertex of the initial 3-clique, its label is first composed of the only digit not 
appearing as first digit of its neighbors (in this case, 1), concatenated with 
the longest label of its neighbors (in this case, 23). Analogously, in the lower 
part of Figure El, where, for sake of clarity, only a part of A(2,3) is drawn, 
the vertices have been labeled. For the same reasons, we can see that vertex 
u has label L(u) = 4, while w has label L(w) = 234. 

Thus, we see that for any t > 1, any vertex vt created at step t has a 
unique label, and that for any vertex v created at step t > 1, L(v) = S1S2 ■ ■ ■ s t 
is of length t, where each digit Sj satisfies 1 < Sj < d + 2 ; while the vertices 
created at step have length 1 (these are the I', 1 < I < d + 2). 

We note that since for any step t > 1, the number of vertices that are 
added to the expanded Apollonian networks is equal to (d + 2)(d + l)* -1 , 
the labeling we propose is optimal in the sense that each label L(vt) of a 
vertex created at step t is a (d + 2)-ary word of length t. Globally, any 
vertex of A(d, t) is assigned a label of length 0(\og d+2 t) ; since there are 
N t = (d+2) {d+1)t d +d - 1 vertices in A(d,t), we can see that, overall, the labeling 
is optimal as well. 

Next, we give three properties about the above labeling. Property 1 
ensures that our labeling is deterministic. Property 2 is a tool to prove 
Property 3, the latter being important to show that our routing protocol is 
valid and of shortest paths. 

Property 1 In A(d, t), for any (d + 2)-clique induced by vertices W\, w 2 
. . . Wd+2, every integer 1 < i < d + 2 appears exactly once as the first digit 
of the label of a Wj. 
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Figure 2: (Above) Labels of all vertices of A(l, 3). (Below) Labels of a part 
of the vertices of A{2, 3) 
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Proof. By induction on t. When t — 1, the property is true by construc- 
tion. Suppose now that the property is true for any t' < t, and let us then 
show it is true for t. Any (d + 2)-clique in A(d,t) is composed of exactly 
one vertex v created at a given step t±, and d + 1 vertices wi,W2, ■ ■ -w^+i 
created at steps strictly less than t\. If t\ < t, then the property is true 
by induction hypothesis. If t-y = t, we suppose that Wd+2 is connected to a 
(d + l)-clique C composed of W\, W2, ■ ■ ■ w^+i- It is clear that C did not exist 
at step t — 1. In other words, one of the w^s, say w±, has been created at step 
t — 1, based on d + 1 vertices W2, W3 . . . Wd+i and x. By induction hypothesis, 
each integer 1 < % < d + 2 appears exactly once as first digit of the labels of 
Wi, w 2 , w 3 . . .Wd + i,x. However, by construction, the first digit of L(u^ +2 ) is 
the first digit of L(x). Thus we conclude that each integer 1 < % < d + 2 also 
appears exactly once as first digit of the labels of w\,W2,w^ . . .Wd+i,Wd+2, 
and the result is proved by induction. 

Property 2 Let v t be a vertex of A(d,t) created at step t > 1. Among 
the vertices w±, w 2 , ■ ■ ■ Wd+i forming the (d + l)-clique that generated v t , let 
Wi,w 2 - ■ -Wk, k < d + 1, be the vertices that do not belong to the initial 
(d + 2)-clique. Then L(v t ) is a superstring of L{wi) for all 1 < % < k. 

Proof. By induction on t. When t — 1, any vertex v\ created at step 1 
is connected to vertices of the initial (d + 2)-clique only. Thus the result is 
true. Now suppose the result is true for any I < t' < t — 1, t > 2, and let us 
prove it is then true for t. For this, we consider a vertex v t created at step 
t, and the (d + l)-clique C it is connected to. Suppose v t is a neighbor of w p 
which was created at step t — 1. However, w p was created itself thanks to a 
(d + l)-clique, say C, composed of vertices x±,X2 ■ ■ -Xd+i- W.l.o.g., suppose 
that k < d+1 such vertices, x\, x 2 ■ ■ ■ Xk do not belong to the initial (d + 2)- 
clique. By induction hypothesis, L(xi) C L(w p ) for any 1 < % < k. Hence, 
in C, w p is the vertex not belonging to the initial (d + 2)-clique that has the 
longest label. By construction of L{y t ), we have that L{w p ) C L{y t ), thus 
we also conclude that L(xi) C L(v t ) for any 1 < i < k. Thus L(v t ) is a 
superstring of the labels of any vertex of C that does not belong to the initial 
(d + 2)-clique, and the result is proved by induction. 

Property 3 Let v t be a vertex of A(d, t) created at step t>l. For any 
1 < i < d + 2, \1% L(vt), then vt is a neighbor of a vertex v' of the initial 
(d + 2)-clique, such that L(v ') = i'. 

Proof. By induction on t. When t — 1, any vertex v 1 constructed at step 
1 is assigned label i, where %' is the only vertex of the initial [d + 2)-clique v t 
is not connected to ; thus, by construction, the property is satisfied. 
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Now we suppose that the property is true for any 1 < t' < t — 1, t > 2, 
and we will show it then holds for t as well. As for the previous property, we 
consider a vertex v t created at step t, and the (d + l)-clique C it is connected 
to. 

Suppose v t is connected to a vertex w t -i that was created at step t — 1. 
However, w t _i was created itself thanks to a (d + l)-clique C composed of 
vertices Xi,x 2 ■ ■ -Xd+i- Among those d + 1 vertices, only one, say x p , does 
not belong to C. W.l.o.g., suppose that k < d + 1 such vertices, x\,x 2 ■ ■ ■ Xk 
do not belong to the initial (d+ 2)-clique. Now suppose that i ^ L(v t ) ; then 
% appears as the first digit of one of the L(xj)'s, j G [1, p — 1] [j[p + 1, d + 1], 
or of L(w t ^i) (by Property 1). However, L(xj) C L(w t _i) C L(v t ) for any 
1 < j < ^ (by Property 2). Thus, neither w t -\ nor any vertex among the 
x/s, 1 < j < k contains the digit i in its label. Hence, only a vertex y from 
the initial (d + 2)-clique can have i in its label, and thus L(y) = i'. Hence it 
suffices to show that v t and y are neighbors to prove the property. The only 
case for which this would not happen is when y = x p ; we will show that this 
is not possible. Indeed, by construction of the labels, the first digit of L(v t ) is 
the only integer not appearing as first digit of the labels of the vertices of C, 
that is Wt-i, xi, X2 ■ ■ ■ x p -i, x p+ i . . . Xd+\- However, the fact that we suppose 
y = x p means that no vertex of C contains i in its label. Thus this would 
mean that the first digit of L{v t ) is i, a contradiction. Thus, v t is connected 
to y with L{y) = i', and the induction is proved. 

4 Routing by shortest path 

Now we describe the routing protocol between any two vertices u and v of 
A(d, t), with labels respectively equal to L(u) and L(v). We note that since 
A(d, 0) is isomorphic to the complete graph K d+2 , we can assume t > 1. The 
routing protocol is special here in the sense that the routing is done both 
from u and v, until they reach a common vertex. Hence, the routing strategy 
will be used simultaneously from u and from v. In order to find a shortest 
path between any two vertices u and v, the routing protocol is as follows. 
First we compute the longest common suffix LCS(L(u), L(v)) of L(u) and 
L(v), then we distinguish two cases: 

1. If LCS(L{u),L{v)) =0: 

(a) Simultaneously from u and v (say, from u) : let u = u and go from 
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Ui to u i+i , i > where u i+ i is the neighbor of u { with shortest 
label. 

(b) Stop when Uk is a neighbor of the initial (d + 2)-clique. 

Let L{uk) (resp. L(vk>)) be the integers not present in L(uk) (resp. 
L(vk')), and let S* = L{uk) fl L{yk>)- 

i. It S ^ 0, pick any / G 5, and close the path by taking the 
edge from Uk to /', and the edge from /' to vy- 

ii. If S — 0, route from -u^ to any neighbor l[ (belonging to 
the initial (d + 2)-clique) of Uk, and do similarly from vy to 
a neighbor 1' 2 (belonging to the initial (d + 2)-clique) of vy ■ 
Then, take the edge from l[ to 1' 2 and thus close the path from 
u to v. 



2. If LCS(L(u), L(v)) 7^ 0, then let us call least common clique of u 
and v, or LCC(u,v), the (d + 2)-clique composed of the vertex with 
label LCS(L(u), L(v)) and the d+ 1 vertices forming the (d+ l)-clique 
that generated the vertex of label LCS(L(u), L(v)). We simultaneously 
route from u and v to (respectively) -u^ and v k r, going each time to the 
neighbor with LCS(L(u), L(v)) as label suffix, and having the shortest 
label. Similarly as above, we stop at u k (resp. tv), where Uk (resp. tv) 
is the first of the u^s (resp. of the v/s) to be a neighbor of LCC(u, v). 
Then there are two subcases, depending on Q = L{uk) H L(vk')- 

(a) If Q ^ 0, close the path by going to any vertex u> with label /, 
leQ. 

(b) If Q = 0, then route from w fc (resp. tv) to any neighbor w 1 (resp. 
W2) in LCC(u, v), and close the path by taking the edge (wi, W2), 
which exists since both vertices w\ and W2 belong to the same 
clique LCC(u, v). 

Proposition 1 The above mentioned routing algorithm is valid, and of 
shortest paths. 

Proof. Let us first give the main ideas for the validity of the above 
routing protocol. Take any two vertices u and v. By construction of L(u) 
and L{v), the longest common suffix LCS(L(u), L(v)) indicates to which 
(d + 2)-clique u and v have to go. We can consider this as a way for u and v 
to reach their least common ancestor in the graph of cliques induced by the 
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construction of A(d,t), or the "least common clique" . In Case (i), this least 
common clique is the initial (d + 2)-clique ; thus, u and v have to get back 
to it. In Case (ii), the shortest path does not go through the initial (d + 2)- 
clique, and the least common clique of u and v, say LCC(u,v), is indicated 
by the longest common suffix LCS(L(u), L(v)). In other words, the length 
of LCS(L(u), L(v)) indicates the depth of LCC(u, v) in the graph of cliques 
induced by the construction of A(d,t). In that case, the routing is similar 
as in Case (i), except that the initial (d + 2)-clique has to be replaced by 
the clique LCC(u,v). Hence, the idea is to adopt the same kind of routing, 
considering only neighbors which also have LCS(L(u), L(v)) as suffix in their 
labels. 

When this least common ancestor is determined, one can see, still by 
construction, that the shortest route to reach this clique (either from u or 
v) is to go to the neighbor which has smallest label, since the length of the 
label indicates at which step the vertex was created. Indeed, the earlier the 
neighbor w was created, the smaller the distance from w to the least common 
clique is. 

After we have reached, from u (resp. from v), a vertex Uk (resp. vy) 
that is a neighbor of the least common clique, the last thing we need to 
know is whether Uk and vy are neighbors. Thanks to Property 3, we know 
that looking at L{uk) and L{vk>) is sufficient to answer this question. More 
precisely: 

• In Case (i)(b)-l, Uk and vy share a neighbor in the initial (d + 2)-clique 
(by Property 3). All those common neighbors have label I', where / G S. 
Hence, if we pick any / e S, then there exists and edge between Uk and 
I', as well as an edge between I' and vy. 

• In Case (i)(b)-2, Uk and vy do not share a neighbor in the initial (d+2)- 
clique. Hence, taking a route from w fc (resp. v k ') to any neighbor l[ 
(resp. 1' 2 ) belonging to the initial (d + 2)-clique, we can finally take the 
edge from l[ to l 2 (which are neighbors, since they both belong to the 
initial (d + 2)-clique) in order to close the path from u to v. 

• In Case (ii) (a) , w fc and vy share a neighbor in LCC(u,v). Hence we 
can close the path by going to any vertex w with label I, I e Q, since 
w is a neighbor of both Uk and tv- 

• In Case (ii)(b), Uk and vy do not share a neighbor in LCC(u, v). Hence 
we route from w fc (resp. v k >) to any neighbor w\ (resp. w 2 ) in LCC(u, v), 
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and we close the path by taking the edge (wi, w 2 ) ■ This edge exists since 
both vertices W\ and w 2 belong to the same clique LCC(u,v). 

Hence we conclude that our labeling of vertices in A(d, t) allows a routing 
between any two vertices u and v, and that it is of shortest paths. 

5 Conclusion 

We have proposed an expanded deterministic Apollonian network model, 
which represents a transition for degree distribution between exponential 
and power law distributions. Our model successfully reproduces some re- 
markable characteristics in many nature and man-made networks. We have 
also introduced a vertex labeling for these networks. The length of the label 
is optimal. Using the vertex labels it is possible to find in an efficient way 
a shortest path between any pair of vertices. Nowadays, efficient handling 
and delivery in communication networks (e.g. the Internet) has become one 
important practical issues, and it is directly related to the problem of finding 
shortest paths between any two vertices. Our results, therefore, can be useful 
when describing new communication protocols for complex communication 
systems. 
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